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ABSTRACT 

The interaction of a nonstationary plasma flow with a model ionosphere is studied. On the 
basis of a numerical simulation, the calculation yields results of the distribution of the 
plasma concentration and magnetic field in the transition region, 

INTRODUCTION 

Three models of the solar-wind interaction with the upper atmosphere of planets not having 
their own magnetic field sufficient for deceleration of the solar wind are known: 

1 . Solar plasma deceleration in the atmosphere, 

2. The boundary between the solar wind and atmosphere is considered as a tangen- 
tial discontinuity, and 

3. The induced magnetosphere excitation [1,2]. 

The first and the second models were studied in detail [3, 4, 5] , but there are only qualita- 
tive estimates of the model with magnetic field induction [6, 7] . 

This paper is an attempt to study the structure of the magnetic field excited by currents in 
the solar wind and the ionosphere. 

FORMULATION OF THE PROBLEM 
The Solar-Wind Model 

The solar-wind parameters are such that >$> D/Dt £n H, v. « W <5< v^ , « p. < 

where p.^ is Larmor radius of ions and electrons, is the radius of the ionosphere, is 
the gyrofrequency of electrons, v.^ and W are the thermal and directed velocity of ions and 
electrons, and H is the modulus of the magnetic field. To describe the electrons, the drift 
approximation is applied, and the electron distribution is thus approximated by a Maxwell- 
Boltzmann function: 
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In the approximation considered, the ions are assumed to be cold and their motion, in 
dimensionless form, is defined by the equation: 

V* = ^ (E* + V* X H*) 
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where and are the plasma velocity and the magnetic field in the incoming flow. 
Maxwell’s equations, in dimensionless form, are as follows [8] : 
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where 

n^ = the plasma concentration in the incoming flow (n* = n/n^ ), 
j3g = /H^ , the ratio of thermal pressure to magnetic pressure, 

X = P||g/Pj^^, the anisotropy of the pressure, 

-v/47rn^ M/H„ , Alfven Mach number, 

K = the curvature of the magnetic field. 

To calculate the electric field, the condition of quasineutrality of plasma n^ <=« n. is used. 
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The Ionosphere Model 

Two models of a planetary ionosphere are considered. In the first model, the ionosphere is 
simulated by a motionless spherical ring < r < r^ with an effective homogeneous iso- 
tropic conductivity a. (See figures 1 and 2.) In this case, the equations of the field are as 
follows: 


V X iT* = R E*, 

e ’ 
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— H*+VX E* = 0 
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R = = the magnetic Reynolds number . 

The boundary conditions on the outer and inner surface of the ionosphere are 

[HJ=[E,]=0, H„=E^=0. 

In the second model, the ionosphere is simulated by a motionless, magnetized, ideally-con- 
ducting, plasma-like, hemispherical cap with a right circular cylinder on the nightside 
(figures 3 and 4). In this model, the characteristics of the ionosphere are defined by the 
parameters 

8nPj. V47m.M 


For simplification, the index * of the dimensionless values will be omitted. 

Initial and Boundary Conditions 

Initially, the solar-wind plasma density and the electromagnetic field in the calculation 
region are absent (in the first model, it is a sphere; in the second model, it is a cylinder). 
Then, the plasma with a frozen-in magnetic field is ejected against the ionosphere. 

On the outer boundary of the calculation region, the parameters of the field and plasma are 
assumed to be undisturbed. In the case where 1 H^, mirror symmetry relative to 
equatorial and meridional planes takes place. Let us expand the parameters in Fourier 
series of 4>, relative to the Z-axis, and consider the homogeneous solution described by the 
first harmonic. 

The solution of equations 1 through 6 was determined by the “guiding center in the cell” 
method, a description of which can be found in Lipatov [8] . 

The solution of Maxwell’s equations was determined by the alternating-direction method 
[9] on the spherical (first model) and cylindrical (second model) grid. 
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Calculation Results 

The calculations were made for the following values of the solar-wind parameters: 

® X = 1 to 1.5, 

® M^=0.3tol, 

® W„/(oj,R^) = 0.01 to 0.1. 

Figures 1 and 2 show the distribution of the normalized plasma concentration n/n^ and 
magnetic field H/H^ near the ionosphere in the meridional and equatorial planes in the case 
of a spherical model of the ionosphere (first model). The directed velocity of the solar-wind 
flow near the nose of the ionosphere reaches ~W^/10, evidence of the plasma flow 
deceleration. The conductivity of the ionosphere is considered to be rather high (R^ ~ 10^) 
so that the time of penetration of the interplanetary field into the ionosphere is high com- 
pared with the time of the particles’ transit past the planet. 

Figures 3 and 4 give the distribution of the normalized plasma concentration n/n^ and mag- 
netic field H/H„ near the ionosphere in the meridional and equatorial planes in the case of 
the existence of the ionospheric tail (second model). The parameters of the ionosphere 
model have the values = 10, = 10^. Drift currents in this model create an intense mag- 

netic field between the solar wind and the ionosphere and screen the ionosphere from pene- 
tration by the interplanetary field. 

In the models of flow considered, the amplification of the magnetic field at the nose of the 
ionosphere is defined, in general, by compression of the lines of force of the field on the 
equator by analogy with plasma flow around the Moon [8] . 

The induced magnetosphere size is defined by the balance of the total pressure of the solar 
wind + 2 p and of the pressure of the induced field However, the 

calculated model of the induced magnetosphere does not take into account all the processes 
of the real deviated flow (see [10]), that is, the situation when the solar- wind ions can 
penetrate the ionosphere. In this case, an additional current system is set up coupling the 
solar wind into the ionosphere. These currents, according to estimates [7], may be impor- 
tant for determination of the magnetic field. 

Filling up the induced magnetosphere by plasma will have an influence on the characteristics 
of the magnetosphere (duration of change of sign of the magnetic field and so on). 

CONCLUSION 

Preliminary calculations of decelerating plasma flow on the model ionosphere show that the 
induced magnetic field, in general, is determined by the conductivity of the ionosphere and 
can reach high values. 
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However, the model of an induced magnetosphere considered does not take into accoimt the 
effect of the closing of currents from the solar wind into the ionosphere and the motion of 
ions of the ionosphere, both of which may influence the interaction of the solar wind with 
the ionosphere. 
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QUESTIONS 

Lipatov iNess: Your diagram shows an increased field near the equatorial plane and a de- 
creased field near the polar region assuming that the equator is defined by the plane con- 
taining the field lines and the solar wind velocity. Can you indicate the geometry of the field 
lines around the obstacle to the plasma flow? Do you expect the increased field to persist as 
the Mach number increases from subsonic to supersonic? 

Lipatov. In the calculated model, we give the configuration of the magnetic lines of force, 
which is similar to the configuration of the lines of force near the Moon. As the Mach num- 
ber increases to a supersonic value, the increased magnetic field remains near the, nose of 
the ionosphere. But in the polar regions we can see an increase of the tangential component 
of the magnetic field. 
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